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Abstract 

We investigate the asymptotic behaviour of spin-spin correlation functions 
for the integrable Heisenberg chain. Tja this end we use the Quantum Transfer 
Matrix (QTM) technique developed inlil which results in a set of non-linear inte- 
gral equations (NLIE). In the case of the largest eigenvalue the solution to these 
equations yields the free energy and by modifications of the paths of integration 
the next-leading eigenvalues and hence the correlation lengths are obtained. At 
finite field h > and sufficiently high temperature T the next-leading eigenvalue 
is unique and given by a I-string solution to the QTM taking real and negative 
values thus resulting into exponentially decaying correlations with antiferromag- 
netic oscillations. At sufficiently low temperatures a different behaviour sets in 
where the next-leading eigenvalues of the QTM are given by a complex conjugate 
pair of eigenvalues resulting into incommensurate oscillations. 

The above scenario is the result of analytical and numerical investigations 
of the QTM establishing a well defined crossover temperature T c (h) at which 
the I-string eigenvalue to the QTM gets degenerate with the 2-string solution. 
Among other things we find a simple particle-hole picture for the excitations of 
the QTM allowing for a description by the dressed charge formulation of CFT. 
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1 Introduction 



In this contribution we report on a new crossover phenomenon observed in the longi- 
tudinal correlation function at finite magnetic field h and finite temperature T of the 
XXZ spin-chain with anisotropy parameter < A < 1. 

As is well known any thermodynamic quantity derived from the free energy of the 
one-dimensional XXZ model is an analytic function of finite h and T . Phase transi- 
tions and associated mathematical singularities may and do occur in the ground state, 
i.e. at T = 0. However, quantities not obtained from the free energy, but characteris- 
ing the asymptotics of correlation functions may show "their own" non-analyticities at 
finite temperatures. Indeed the XXZ chain in an external magnetic field does show 
well defined non-analyticities in the correlation lengths of the longitudinal spin-spin 
correlation functions. We want to point out that these crossover phenomena are a 
result of strong correlations and finite temperature. The crossover does not take place 
for the case A = corresponding to free fermions. 

For the attractive regime — 1 < A < of the XXZ chain the investigations are 
not yet carried out for finite field h and finite T. Here the physical phenomena are 
expected to be much richer. Even for the case of vanishing magnetic field h — a se- 
quence of crossovers (commensurate-incommensurate-commensurate) for an increase 
of T from to oo was found □. n 

The method we apply (as that used incl) is based on a mapping of the XXZ chain 
at finite temperatures onto a classical model in two dimensions and the formulation 
of a suitable quantum transfer matrix (QTM) describing the transfer along the chain. 
The spectrum of the QTM succumbs to a Bethe ansatz (BA) treatment and the cor- 
responding BA equations can be cast into the form of non-linear integral equations. 
Crossover phenomena of correlation lengths manifest themselves as level crossings 
of next-leading eigenvalues of the QTM. At finite field h > and sufficiently high 
temperature the relevant eigenvalues for the longitudinal correlation functions are 1- 
string and 2-string solutions (both solutions belong to the S = sector of the model) . 
The truely next-leading eigenvalue is unique and given by the f-string solution to 
the QTM taking real and negative values thus resulting into exponentially decaying 
correlations with antiferromagnetic oscillations. In some sense at sufficiently high 
temperature the properties of the system are determined by the longitudinal ("clas- 
sical") terms of the Hamiltonian, i.e. the S Z S Z coupling and the field in z direction, 
which dominate over the transversal exchange ("quantum mechanical") terms. At suf- 
ficiently low temperature a different behavior is expected on grounds of predictions 
by conformal field theory (CETpQ. In particular, correlations with incommensurate 
2k f oscillations are expectecErtl As a consequence of this, the QTM has to develop 
complex conjugate pairs of eigenvalues at sufficiently low temperatures. This scenario 
has not been investigated before. The purpose of this contribution is to present the 
first "Bethe ansatz" study of this crossover phenomenon and to provide reasonably 
accurate values for the crossover temperature T c . As pointed out already, for the free 
fermion case T c = oo, i.e. the crossover does not take place. In physical terms this 
may be understood in the way that due to the absence of longitudinal couplings the 
longitudinal terms never dominate over the transversal terms. 

This report is organized as follows. In Sec. 2 we introduce some basic definitions 
of the XXZ chain and present its properties at low temperature as obtained within 
conformal field theory (CFT). In Sec. 3 the approach to thermodynamic properties 
by use of a lattice path integral formulation and the quantum transfer matrix (QTM) 
is reviewed. In particular the set of non-linear integral equations for the two auxiliary 
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functions a(x) and o(x) corresponding to the energy density functions of spinons with 
spin ±1/2 are given. In Sec. 4 numerical results are given for the correlation length 
and Fermi momentum of the longitudinal spin-spin correlation function. Finally, the 
particle-hole picture resulting at low temperatures is discussed. A complete exposition 
of how this is related to the dressed charge formulation of CFT will appear int3. 



2 Anisotropic Heisenberg model 

The XX Z model is defined by the Hamiltonian 

// -I £ [S?Sf + SfS* + ASfS?] h^Sj, 

<i,j> j 

where J > 0, h > and S x ' v,z are spin-1/2 operators. For T = the system is in one 
of three phases. For anisotropy parameter |A| > 1 we have two ordered phases and for 
|A| < 1 a critical phase. In the latter case the anisotropy parameter is conveniently 
parameterized by A — cos 7. In the following we take < 7 < tt/2 (repulsive regime) 
for simplicity. The excitations are gapless spinons with spin 1/2 

e F (k) = v smk, < k < tt, (2) 



v = nJ. (3) 

7 

At zero temperature the spin correlations decay algebr^call^S00lll, with expo- 
nents depending on the magnetic field in a nontrivial wayl3. At finite temperature the 
correlations decay exponentially e~ r /£ with in general different correlation lengths £ 
for different correlation functions. At low temperatures the correlation lengths can 
be related by conformal mappings to the scaling dimensions x of the fields at T = 

£ = for T«l 

5 2irxT 

The scaling dimensions x are calculated from finite size corrections to the energy 
levels of the Hamiltonian and scaling predictions by CFT 

E x -E = ^v(x + N+ + N-) + o^y (4) 

where N + and N~ are integers labelling the conformal tower. The critical exponents 
of the spin-1/2 XX Z model are those of a c = 1 Gaussian theory. At zero magnetic 
field the exponents are given in terms of the anisotropy parameter 7 

where S and m are integers corresponding to the S z component of the state and 
the number of excitations from the left to the right Fermi point, respectively. The 
longitudinal correlation is 

,ozoz\ ,, n cos(2fc r r) C 2 , , 

\ b o b r) ~ 0l r i/(i- 7 / w ) " 72' W 

where C\ and C2 do not depend on the distance r. For zero magnetic field the Fermi 
momentum k F is equal to 7r/2. For finite magnetic field and temperature we expect 
a deviation of x and kp from the above quoted values. 
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3 Finite Temperatures 



The properties of the quantum system at finite temperature are determined within a 
path integral formulation. The partition function of the Id quantum system at finite 
temperature is mapped to that of a 2d classical system pictorially represented by 



, QTM 



Tre -/3ff =par t.fct. 




time 



For the classical system we can choose a suitable transfer approach. The quan- 
tum transfer matrix (QTM), i.e. the column-to-column transfer matrix, possesses 
a unique largest eigenvalue and a spectral gap to the next-leading eigenvalues which 
persist s eve n in t he limit of infinite Trotter number as long as the temperature is 
finitJaOraBEEl . Hence the free -energy and correlation lengths are simply deter- 
mined from the largest Aq and next largest eigenvalues Aj of the QTM 
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with coefficients A%, A2 given by matrix elements. The diagonalization of the QTM 
can be performed by means of a Bethe ansatz resulting into eigenvalue equations of 
the form of Baxter's A — q equation. For finite Trotter number N we find 

A(x)q(x) = e~ ph/2 ^(x - i"f/2)q{x + ij) + c +(ih/2 <5>{x + n/2)q(x - ij) (8) 



where 



<f>(.x) = [sinh(x — ixo) sinh(a; + ixo)] N ^ 2 , 
q(x) = J^smh(x — Xj) 



1 a 



(9) 



The roots determined from the Bethe ansatz equation 



p{xj) = — 1, where p(x) := c 



$(x + ij/2)q(x — ij) 



(10) 



thus rendering A(a;) analytic. We like to note that in general there are more solutions 
to the Bethe ansatz equation p(x) = —1 than roots Xj. The additional solutions are 
called holes. For an illustration of a typical case see Fig [j]. 

The corresponding distributions of roots and holes remain discrete even in the 
limit of infinite Trotter number N — > 00. The method of studying these equations 
was developed inu and results into two equations for two auxiliary functions a(x) := 
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Figure 1: Distribution of roots and holes for the largest eigenvalue for anisotropy parameter 7 = 7r/3, 
reciprocal temperature /3 = 12 and h = 0.6. Note that straight contours with imaginary parts ±7/2 
separate the roots from the holes. The contours C\ — if/2 and £2 + if/2 are depicted by dashcd- 
dotted lines. Note that roots and holes get arbitrarily close to C\ — if/2 for low temperatures T and 

positive h. 



l/p(x— 27/2) and a(x) := p(x+i~//2) (corresponding to S = 1/2 spinon and antispinon 
dressed energy functions) 

lno(x) = h n , 7r ^ h , + [k *i ln(l + a)](x) - [k * 2 hi(l + a)] (a: - 17 + ie). 

cosn —x 2[tt — 7) 

(11) 

An analogous integral equation for a can be obtained from the above one by use of 
the obvious identity lna(ir) = — lna(ir + 17) thus completing the non- linear integral 
equations. The integral kernel k(x) is defined by a Fourier integral 

k(x) := i- f Sinh(f ~ 7)fcC0S(M dfc. (12) 
v ; 2cosh2jfc sinh^/c v ; 

The symbol *c denotes convolution 

/ *g{x) = J f{x-y)g{y)dy (13) 

with a suitably defined integration contour L. In ( |ll|) the subscripts *i (* 2 ) refer 
to integration paths C\ and £2 that extend from —00 to +00 and lie below the 
distribution of numbers xj + ij/2 and above Xj — *7/2. An additional and last 
requirement for these paths is that all roots Xj be situated between C\ — ij/2 and 
£2 + *7/2, but no hole solution to the Bethe ansatz equation. 

For the largest eigenvalue the distribution of roots and holes is depicted in Fig. [j]. 
Here the paths C\ — ij/2 and C2 + ij/2 are just straight lines with imaginary parts 
T7/2- Hence, £1,2 coincide with the real axis. 
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Figure 2: Distribution of roots and holes for a 1-string solution for anisotropy parameter 7 = 7r/3, 
reciprocal temperature = 3.0, and (a) h = and (b) h = 0.6. Of all hole solutions only the two 
new holes close to the real axis are shown. Note that for (a) the 1-string is situated exactly at — iiz/2 
and the holes lie on the real axis. For positive magnetic field (b) the 1-string is shifted upwards (see 
also Fig. W). The contour L\ — is depicted by a dashed-dotted line. 
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Figure 3: Distribution of roots and holes for a 2-string solution for anisotropy parameter 7 = 7r/3, 
reciprocal temperature /3 = 3.0, and (a) h = and (b) h = 0.6. Only the two new holes close to the 
real axis are shown. Note that for (a) the 2-string is symmetric with respect to the real axis and 
situated approximately at ±«7/2; the holes lie on the real axis. For positive magnetic field (b) the 
2-string is shifted downwards (see also Fig. W). The contour C\ — ry/2 is depicted by a dashed-dotted 

line. 
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The cases of next-leading excitations which are of 1-string type (root yo in the 
lower half plane) and 2-string type (upper and lower roots y+ and y_ separated by 
approximately ij) are shown in Figs. g, [|. For the excited states the contours of 
integration are deformed. Only the path L\ — i~f/2 is shown explicitly, C2 + ij/2 is 
simpler as it is mostly following a straight line with imaginary part +7/2 with loop 
in counter clockwise manner around yo + *7 and y~ + i"f, respectively. 

The contours may be straightened and by use of Cauchy's theorem we pick up 
residues that lead ±0 additive contributions to the driving terms in the non-linear 
integral equations □. This form of the NLIE is particularly useful for numerical 
treatments!]. 

Finally, the subsidiary conditions a(yo + £7/2) = —1 and a(y± + i-f/2) = — 1 
with yo, y+ and y_ denoting the root of the 1-string and the upper and lower con- 
stituents of the 2-string, respectively, yield the information on the positions of the 
string parameters. 

The eigenvalues of the QTM are expressed in terms of a(x), a(x) 

1 f]n[(l + a(x))(l + a(x))] J 
= +2-J coshL — (W) 

where eo is the groundstate energy for zero magnetic field and the integration contours 
are £1,2- However, due to certain analyticity properties of the integrand these contours 
can be simplified to just one contour along the real axis, surrounding the numbers 
8 + ij/2 in clockwise manner where 8 is any of the hole solutions to the Bethe ansatz 
equation close to the real axis. 

The range of validity (convergence) of the integral equation (|lTJ) for a(x) is the 
strip Im(x) G [0,7]. Sometimes it is necessary to extend the above equation to the 
strip Im(x) £ [— (n — 7),0], e.g. for calculating a(yo + *7/2) = — 1 with y in the lower 
half plane, see Figs. One version of such an expression for the 1-string pattern is 

lna(x) = — [r * ln(l + 0)] (x) + [r * ln(l + 0)] (x — i*y) 

7T ~ 7 



log 



sinh^(x-yo-»3 7 /2) sinh ^(x - 8 3 + » 7 /2) 
sinh ^(x -y Q + if/2) .1^ sinh ^(x - 8 3 - i 7 /2) 



(15) 



where 



r ( x ) — ~i 7 coth x — coth (x + 17) . (16) 

2(TT — 7) \ 7T — 7 7T — 7 J 

Note that the above expression for a(x) is valid for x below the real axis, however 
a, a on the right hand side are evaluated for strictly real arguments, i.e. * denotes 
standard convolution with the integration contour along the real axis. In particular, 
there are no deformations of the contour as discussed above. They have been removed 
by use of Cauchy's theorem. In turn, the rapidities 0±, 82 and the 1-string parameter 
yo show up explicitly. Note that the corresponding expression for the 2-string pattern 
is obtained by simply exchanging yo with y_ (the lower part of the 2-string). In this 
case, the upper part y + does not appear explicitly. The details will be presented inl3. 

4 Crossover scenario and the particle-hole picture 

From the next-largest eigenvalue(s) Ai of the QTM the asymptotic behaviour of the 
correlation functions is determined. If Ai is real and positive (negative), the decay 
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Figure 4: Schematic depiction of the crossover scenario. At the left we show for the 1-string and 
2-string solutions how the strings move away from their high temperature positions upon a decrease 
of T. At some well-defined temperature T c both solutions become degenerate. A further decrease of 
temperature forces the f-string and the lower part of the 2-string to move away from the imaginary 
axis. As a result we have two complex conjugate solutions with complex eigenvalues. 



is purely exponential without (with) sublattice oscillations. In general Ai may be 
complex and the asymptotics are characterized by a finite correlation length £ and 
Fermi-momentum Isf defined by 

Al _ -l/£±i2k F 

A 

We have numerically solved the nonlinear integral equations for the next-leading eigen- 
value^) and plot the results in Figs. ||, ||. 

Summarizing our results we see that the previously drawn picture of excitations of 
distinctly 1-string and 2-string type is valid only for h = or, if h > for sufficiently 
high T. For this case the 1-string solution dominates and the eigenvalue is real 
because of the left-right symmetry of the Bethe ansatz pattern. The 2-string solution 
is subdominant. Both strings lie on the imaginary axis where for T = oo the 1- 
string is located at —itt/2 and the 2-string is symmetric with respect to the real 
axis. For decreasing temperature we observe a characteristic motion of the 1-string 
upwards and the 2-string downwards along the imaginary axis, cf. Fig. ||. This motion 
continues until the root of the 1-string (y ) and the lower root of the 2-string (y_) take 
identical values. The corresponding temperature defines T c . For lower temperatures 
a horizontal motion sets in. The considered root previously on the imaginary axis 
develops a non-vanishing real part which may be positive or negative, see Fig. |]. 
The corresponding Bethe ansatz patterns are related by reflection at the imaginary 
axis, the corresponding eigenvalues are complex conjugate. There is no longer any 
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qualitative distinction like for temperatures higher than T c . 

The reason for this crossover may be understood qualitatively in the following 
way. The subsidiary condition for the 1-string j/o ( or f° r the 2-string with tjq replaced 
by is a(yo + *7/2) = —1 which yields due to ( |f5| ) 



log 



7T — 7 

= "integral expressions". 



sinh -^(yo -9i + ij) sinh -^(y -6 2 + ij) 
sinh^^o-^i) sinh^(y -e 2 ) 



(17) 



Unfortunately, the right hand side has to be evaluated numerically. It turns out to 
be of order 0((3h), but smaller than the first term on the left hand side. In any 
case, the equation ( fTi] ) for yo has two inequivalent solutions. These solutions are 
purely imaginary and different for small (3h (<-> distinct 1- and 2-strings); they have 
same imaginary part but non- vanishing real parts with opposite signs for large (3h (<-> 
degenerate 1- and 2-strings). The crossover is typically associated with square root 
singularities. The precise value of T c can only be calculated numerically 

The results for the correlation length and Fermi momentum are shown in Fig. g. 
Most strikingly we see a non-analytic temperature dependence for non-vanishing ex- 
ternal magnetic fields at the well-defined crossover temperature T c \ The singularity 
at T c is of square root type. Most significant is the non-analytic behaviour of the 
Fermi momentum. At low temperature the oscillations are incommensurate and at 
zero temperature the Fermi momentum kp and magnetization m are strictly related 
by fe_p = (1/2 — m)n, a relation which ceases to hold at elevated temperatures. At 
sufficiently high temperatures the oscillations are commensurate with fcj? = ir/2. The 
loss of the left-right symmetry in the Bethe ansatz patterns at low temperatures is 
the reason for the incommensurability of kp and the non-analytic behavior of the 
correlation length. 
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Figure 5: The XXZ chain for finite field h = 0.1 and several values of the anisotropy parameter 7. 
Temperature dependence of: (a) longitudinal correlation length £ times T, (b) Fermi- momentum. 
The crossover temperature is monotonously decreasing for decreasing 7 with a seemingly finite lim- 
iting value for 7 — * 0. 

We note that temperature and magnetic field act in roughly opposite ways which 
can be inferred in part from their appearance in the combination h/T in the NLIE. At 
sufficiently low temperatures the distinguishing characters of the 1-string and 2-string 
disappear completely and a root-hole picture emerges as shown in Fig El Here the 
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relevant modifications of the excited state of the QTM in comparison to the ground 
state (see Fig. [I]) can be characterized as a rearrangement of the roots and holes on 
an ellipsoidal curve. Clearly, the same constructions and terminology as used for the 
excitations of Fermi systems applies. In this finite field case (h > 0), the analytic 
treatment of the low-temperature asymptotics not only confirms the predictions of 
CFT, but also recovers completely the dressed charge formalism as will be shown irIZl. 
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Figure 6: Bethe ansatz pattern at low temperature (h = 0.6 and (3 = 12.0) corresponding to a 
1-string (or 2-string) solution at sufficiently high temperature. 
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